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ABSTRACT 

We measure the merger fraction of Type 2 radio-loud and radio-quiet active galactic nuclei at z > 1 
using new samples. The objects have HST images taken with WEC3 in the IR channel. These samples 
are compared to the 3CR sample of radio galaxies at z > 1 and to a sample of non-active galaxies. We 
also consider lower redshift radio galaxies with HST observations and previous generation instruments 
(NICMOS and WEPC2). The full sample spans an unprecedented range in both redshift and AGN 
luminosity. We perform statistical tests to determine whether the different samples are differently 
associated with mergers. We find that all (92%+®^^) radio-loud galaxies at z > 1 are associated 
with recent or ongoing merger events. Among the radio-loud population there is no evidence for any 
dependence of the merger fraction on either redshift or AGN power. For the matched radio-quiet 
samples, only 38%1}5 are merging systems. The merger fraction for the sample of non-active galaxies 
at z > 1 is indistinguishable from radio-quiet objects. This is strong evidence that mergers are the 
triggering mechanism for the radio-loud AGN phenomenon and the launching of relativistic jets from 
supermassive black holes. We speculate that major BH-BH mergers play a major role in spinning up 
the central supermassive black holes in these objects. 

Subject headings: galaxies:active — galaxies:interactions — galaxiesrjets — galaxies:nuclei - X- 
rays: galaxies 


1. INTRODUCTION 

One of the most important issues in modern as¬ 
trophysics is understanding the co-evolution of galax- 
ies and their cent r al supermassive black holes (SM BH) 
(jHeckman fc BestI 12014 [Alexander fc HickoxI 120121 for 
recent reviews on the subject). Both numerical sim¬ 
ulations and theoretical arguments show that black 
hole (BH) growth occurs during short-lived periods 
10^ — 10®yr) of intensive accretion whic h are also as¬ 
sociated with powerful Quasar activity dSpltan JlO^ ; 
Reeslll984 iDi Matteo et ani2005l: iHopkins et al.l 2008b : 
Somerville et al.l 1200811 . These are also expected to cor¬ 
respond to periods in which galaxies grow hierarchi¬ 
cally. Since the matter that ultimately accretes onto the 
central black hole needs to lose almost all 99.9%) 
of its angular momentum, studies of mergers, tidal in¬ 
teractions, stellar bars and disk instabilities are cen¬ 
tral for understanding the details of such a process. 
Numerical simulations and analytic calculations have 
shown that major (gas-rich) mergers are capable of 
efficiently driving gas inflows towards the central re- 
gion of the galaxy (< Ikpc) through t i dal forces (e.g . 
Hernauisfl 119891: [Pi Matteo et al.l 120051 iLi et al.l l2007t 
Hopkins fc Quataerfll2011h . and ultimately drive the gas 
to the central ^ Ip c, forming an accretion dis k around 
the central SMBH (jHopkins fc QuataertJl20f(ih . This is 
also thought to be a likely scenario for the formation 
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of galaxy spheroids (e.g. iHopkins et al.l[^08all . How¬ 
ever, disk instabilities and minor mergers may also be 
able to provide material for black hole ac cretion (e.g. 
iHernauist fc Mihoill995t iMenci et al.ll2014ll . 

It is extremely important to study the connection be¬ 
tween galaxy and black hole growth from a purely ob¬ 
servational point of view. One of the central questions 
is whether mergers or other mechanisms may constitute 
the main triggering mechanisms for active galactic nu¬ 
clei (AGN). A number of recent papers investigated this 
issue using data from different surveys of galaxies and 
AGNs. Results are often contradictory. Although it is 
known that not all of them are AGNs, powerful Ultra- 
Luminous Infrared Galaxies ( ULIRG) are ubiquitously 


Veilleux et al. 

1200211. According to some models (e.g 

Hookins et al 

I2008bllal) these objects are believed to rep- 


resent a fundamental stage in the process of the for- 
mation of elliptical galaxies. iKartalteoe et al.l (120121 1 
showed that for a sample of z ^ 2 UL IRGs in the 
GOODS-South held (|Giavalisco et al.ir2004 l. the fraction 
of mergers is up to 70%. For lower luminosities ob¬ 
jects (LIRGs) in the same redshift range the merger frac¬ 
tion found by the same authors is signihcantly smaller 
(- 30%). 

While the merger-ULIRG connection seems to be well 
established, whether this is a vi able scenario fo r all A GNs 
is still an unanswered question. iBahcall et al.l (119971) ob¬ 
served a sample of 20 relatively nearby QSOs with the 
Hubble Space Telescope (HST) and found that the ma¬ 
j ority of them reside in merging galaxies. IGrogin et al.l 
(|2003f) studied the optical counterparts of the X-ray se- 
lected AGNs from th e IMsec Ghandra Deep Field South 
(|Giacconi et al.l[2002h . Based on both the asymmetry in¬ 
dex and the frequency of close companions, these authors 
concluded that mergers and interactions are not good 
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indicators of AGN activity. Interestingly, using HST 
images taken with the Wide Field Camera 3 (WFC3) 
iSchawinski et al.l (|2012f) showed that only a small frac¬ 
tion of h eavily obscured QSO hosts are associated with 
mergers. iTreister et al.l (j2012fl put together a relatively 
large sample of literature data and found clues for a lumi¬ 
nosity dependence of the merger fr action in both ULIRGs 
and AGNs. On the other hand, iVillforth et al.l (1201411 
showed that their sample of QSOs observed with HST 
did not show any evidence for a dependence of the merger 
fraction on lu minosity. Recently, b ased on a sample of 
SDSS galaxies. ISabater et al.l (|2015[ 1 found that the effect 
of interactions is minimal in triggering AGN activity. 

A major issue is related to the so-called radio- 
loud/radio-quiet dichotomy of active nuclei. It has been 
long argued that in order to produce powerful relativis¬ 
tic jets, radio-loud AGNs must possess an extra source 
of energy with respect to radio-quiet AGNs. The most 
popular scenarios among those proposed so far assume 
that the energy may be extracted from the innermost 
region of a magneti zed accretion disk around a rapidly 
spinning black hole iBlandford fc ZnaieH (|1977f) . In the 
Blandford-Znajek framework, the radio-quiet/radio-loud 
dichotomy can be explained in terms of a correspond¬ 
ing lo w/high black hole spin separation (jBlandford et al.l 
llQQOf) . Alternative models pre dict that jets are directly 
powered by the accretion disk (IBlandford fc Pavnd[l982l : 
iPudritz fc NormanI 119861 : iXu et al.l 1199911 . Both theo¬ 
retical arguments and observational e vidence support 
the B landford-Znajek scenario (see e.g. IGhisellini et al.l 
1201411 . Recent numerical simulations also seem to con¬ 
firm such a mechanism as a via ble physical explanation 


for efficient jet production fe.g. iHawlev fc KrolikI 


McKinne^^Gammy 12004 iTchekhovskov et al.l 


2006 


2011 


Sadowski et al.l 1201511 . iWilson fc GolbertI ( 19951 1 origi¬ 


nally proposed that radio-loud AGNs are associated with 
rapidly spinning black holes that are ultimately spun-up 
by major BH-BH mergers. It is clear that a straightfor¬ 
ward prediction of such a scenario is that radio galaxies 
and radio-loud QSOs (RLQSOs) should exhibit signa¬ 
tures of major galaxy mergers. 

Radio-loud AGNs are known to be hosted by large 
elliptical galaxies and are associ a ted w i th SMBH of 
at least ^ 10^ Mr^ (e.g. iLaod 120001: [Punlop et al.l 

l2003l:[Best et al.ir2004 iChiaberge fc Marconill2011[l . Fur- 

thermore, these objects are present in richer envi¬ 
ronments (cluste rs and groups o f galaxies) than ra¬ 
dio quiet AGNs (IShen et al.l 120091 : IDonoso et al.l 120101 : 
iRamos Almeida et al.l 1201311 and they are often associ¬ 
ated with brightest cluster galaxies (|Best et al.l 1200711 . 
At least a fraction of them have been known to be 
associated with merging system for a long time (e.g. 


Heckman et al.lll986HGolina fc de ,luanlll995h . Recently, 
Ramos Almeida et al.l ()2012| 1 studied samples of rela¬ 


tively low-redshift {z < 0.7) sources with deep ground 
based observations. These authors found that the large 
majority (^ 80%) of radio-loud objects show disturbed 
morphologies. The same group also studied a small sam- 
ple of radio-quiet Type 2 quasars (|Bessiere et ahllMT^ 
and found that the merger fraction sample is as high as 


75%. 

However, an accurate census of the merger fraction in 
carefully selected samples of AGNs over a large range 
of redshift, luminosity and radio-loudness is still needed 


to provide a final answer to the above questions. This 
should also be based on a homogeneous set of deep, high 
spatial resolution observations and supported by firm 
statistical evidence for any difference in the observed 
merger fractions amongst the different samples. The aim 
of this paper is to make a significant step towards such 
a goal. 

With the aim of determining the importance of merg¬ 
ers in triggering different types of AGN activity, we se¬ 
lect samples of both radio-loud and radio-quiet Type 2 
AGN. In this paper we focus on Type 2 objects only, be¬ 
cause the bright nuclear component that dominates the 
emission on Type 1 AGNs (QSOs) hampers a detailed 
morphological study of both the host galaxies and the 
close environment of those objects. This constitutes a 
serious concern, particularly at moderate to high red- 
shifts, even using HST images. In this work we focus on 
objects with 1 < z < 2.5. This is a range of redshift 
where the peak of the AGN activity is believed to oc¬ 
cur, and where mergers play a dr amatic role in the lat e 
evolution of massive galaxies (e.g. iHookins et al.ll^OSbll . 
The availability of HST data for a substantial number of 
objects in each sample is key to this work. In particular, 
we focus on samples observed with WFC3 and the IR 
channel. The extremely high sensitivity, low background 
level and the range of wavelengths covered by such an in¬ 
strument is particularly suitable to our goals. In fact, the 
rest frame wavelength corresponding to the most widely 
used WFC3-IR filters (i.e. F140W and F160W) is still 
well within the optical range for objects in the redshift 
range of our interest. Even in relatively short exposures 
(i.e. a few hundreds seconds, or less than 1 HST orbit) 
the effects of merger events on the structure of galaxies 
are clearly revealed in WFG3-IR images. This is possible 
thanks to the high sensitivity in the near IR coupled with 
the large field of view of WFC3, compared to previous 
generation instruments such as WFPC2, NIGMOS and 
AGS. 

In Sect.[2]we describe the samples and the observations 
analyzed in this work, while in Sect. [3] we discuss our 
method to classify mergers based on visual inspection of 
the HST images. In Sect. |4]we describe the statistical 
analysis of the results. In Sect. [S]we discuss our findings 
and we outline our framework for their interpretation in 
Sect. [6] Finally, in Sect. 0 we summarize our results and 
we draw conclusions. 

Throughout the paper we assume = 0.3, = 0.7 

and Hq = 70 Km s“^ Mpc“^. For the magnitude system, 
we use AB magnitudes. 

2. THE SAMPLE 

One fundamental goal of this work is to establish 
whether there is a significant difference in the role played 
by mergers among different classes of AGNs and for dif¬ 
ferent bins of redshift and bolometric power. To this 
aim, we use well defined samples of classical radio galax¬ 
ies spanning 5 dex in radio power. We also select new 
samples of both radio quiet and radio loud AGNs, as well 
as non-active galaxies at z > I. In doing so, we use par¬ 
ticular care in separating the radio quiet and radio loud 
populations, as explained below. 

In the following we describe the properties of each of 
the samples used as part of this work, together with the 
selection criteria used to derive new samples. 
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2 . 1 . Radio loud AGNs: z > 1 SCR radio galaxies 

The first sample we consider is compose d by the radio 
galax ies with z > 1 from the 3CR catalog (jSDinrad et al.l 
1198511 . These are all powerful radio g alaxies belonging 
to th e Fanaroff-Riley class II (FR II iFanaroff fc Rile\i 
[1971 . i.e. those in which the brightest components of 
the radio structure lie at the edges of the radio source, 
in contrast to the FR Is, in which the peak of the radio 
emission is located at the core. 

The 3CR is a flux limited sample selected at low ra¬ 
dio frequencies (Siyg > 9 Jy at 178MHz). Since the 
radio emission at such frequencies is dominated by the 
radio lobes, the selection is independent of the AGN 
orientation. The original catalog includes both Type 1 
(QSOs) and Type 2 objects (radio galaxies) and it is 
one of the best studied samples of radio loud AGNs, be¬ 
ing perfectly suitable to test AGN unification scenarios. 
The high redshift objects included in the 3GR are all 
firmly established radio loud AGNs with powerful rel¬ 
ativistic jets emanating from the central supermassive 
black hole. Even at the lowest lumi nosities, the active 
nucleus is strongly rad io loud (see e.g. lSikora et al.ll2n07t 
iChiaberge et al.ll2005ll . a,ssuming the canonical th reshold 
R = F^ghz/Fb > 10 (iKellermiann et al1ll989tl . where 
F^GHz and Fb are the fluxes in the radio band at 5 GHz 
and in the optical B band, respectively. The z > 1 sub¬ 
sample of the 3GR catalog includes 58 objects. The high¬ 
est redshift object is the radio galaxy 3C257 at z = 2.47. 

The radio power of these AGNs is L 151 ^ 10^® erg 
Hz“^ or slightly higher, which corresponds to a 
bolometric luminosity Lf,oi ^ iq 45-46 gj.g g-i assum¬ 
ing standard bolometric corrections. This result is also 
supported by the estimates of the X-ray luminosity for 
some of these objects, w hich is typica l ly in the range 
-^ 2-10 10 ^ 4-45 gj-g g-i iSalvati et al.l (1200811 observed 

a sample of the most luminous and most distant radio 
galaxies and QSOs from the 3CR catalog with XMM. 
They found that the intrinsic X-ray luminosity is in the 
range 6 x lO^"' < L 2 -iokeV < 2 x 10"^® erg s“^, the Type I 
QSOs being a factor of ^ 6 or higher brighter than the 
(Type 2) radio galaxies. These authors interpreted such 
a discrepancy as a result of the p resence of a beamed 
component in the Type I QSOs. iWilkes et al.l (|20I3[1 
found similar results based on Chandra data. Torres! et 
al. (priv. comm.) analyzed both Chandra and XMM 
archival data for the 3CR radio galaxies observed with 
HST, and found intrinsic X-ray luminosities as low as 
L 2 -wkeV ~ 2 X 10"^^ erg s“^. However, the absorbing 
column density is poorly constrained, and it is possible 
that at least some of the Type 2 radio galaxies are in fact 
Compton-thick. In that case, the derived X-ray luminos¬ 
ity sets a lower limit to the intrinsic AGN luminosity in 
that band. 

Twelve 3CR radio galaxies were observed with WFC3- 
IR and the FI40W filter (in addition to WFC3-UVIS 
and F606W, not used in this work) as part of program 
SNAPI3023 (PI M. Chiaberge, Hilbert et al., in prepa¬ 
ration). The wide-band filter used for the HST obser¬ 
vations (FI40W) includes emission lines at the redshift 
considered here (mainly Ha, H/3 and/or [OHI]5007, de¬ 
pending on the redshift of the source). However, the only 
object that is signi ficantly contamin ated by line emis¬ 
sion is 3C 230 ('see iSteinbrind l2?)TTi l. The HST image 



Fig. 1.— Radio loud samples considered in this paper. The 
horizontal line represents the canonical separation between low and 
high power radio galaxies (i.e. FR I and FR II, respectively). The 
color code identifies the different catalogs, and the shaded area 
show the regions of the radio luminosity v redshift plane where HST 
observations are suitable for this paper (see text). Here, the z< 0.3 
range is covered by the 3CR sources (yellow), and the four groups 
that are part of the Willott sample (6C, 7C and TOOT are in red, 
blue and green, respectively plus the SCR with 2 ~ 0.5) cover the 
range around z=0.5. Bot h the low-z SCR and the Willott sample 
are described in Sect. [TTI The 2 > 1 range is covered again by the 
SCR at high luminosity, and by the ECDFS HzLLRGs (purple) at 
the lowest luminosities. See Sect. [2] for more details. 

of this specific object only shows the narrow-line region 
and thus we cannot observe the stellar component of the 
host. Therefore we exclude this source from our sam- 
pl43. HST snapshot surveys of complete samples are well 
suitable for statistical studies, since the observations are 
scheduled by randomly picking objects from the original 
target list to fill gaps in the HST schedule. The observed 
sample spans the entire range of redshift of the original 
list, i.e. from z = I.O to 2.47. The comparison samples 
described in the following are tailored to match the prop¬ 
erties of these 3CR galaxies, in either bolometric power 
or redshift range (or both). 

In Tab. d] we report the data for the II sources belong- 
ing to this sample. The 24/rm luminosities taken from 
iPodigachoski et ah! (I20I5D are used in the following to 
check that the relevant comparison sample is correctly 
matched to the Hz3C. 

In Fig. [T] we show the location of the 3CR sample in 
the radio power vs. redshift plane, with respect to other 
samples used throughout the paper. Note that, at any 

® A quick inspection of the images from the HST-SNAP13023 
program shows that 3C 230 is the only object that presents an 
identical structure in both of the observed bands which is easily 
recognizable as due to emission line filaments. All other objects 
show, in the filter used in this work, a much smoother morphol¬ 
ogy that is not typical of emission line regions. The bandwidth 
of F140W is extremely wide (~ 0.4/im), therefore the expected 
emission line contamination is small. In fact, for all radio galax¬ 
ies except for 3C 230, we estimate that, based on the observed 
emission line flux from ground based spectra, the total emission 
line contamination does not exceed 20% of the total observed flux 
(Hilbert et al., in preparation). Similar considerations hold for the 
other samples used in this work. 
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redshift, SCR objects are always the most powerful radio 
sources. 

2.2. Radio loud AGNs: low-luminositv radio qalaxies at 
l<z<2.5 

The sample of high redshift low-luminosity radio galax¬ 
ies (HzLLRGs) is derived from the Extended Chandra 
De ep Field South (ECD FS iLehmer et al.ll2005l) . We use 
the IBonzini et al.l (|2012fl catalog and AGN type classifi¬ 
cation to derive this sample. The criteria are as follows. 

i) The spectroscopic, if available, or photometric red¬ 
shift of the source must be in the range 1 < z < 2.5. 
The redshift range is chosen to roughly match the range 
spanned by the SCR sample described above. 

ii) The total radio power at 1.4 GHz must be greater 
than Piaghz > 10^° erg Hz“^ and below the fiducial 
FRI/FRII separation p inHy. < 4 x 10^^ erg Hz“^ 
(|Fanaroff fc Rile\^ll974l l. 

i ii) The source must b e classified as a radio loud AGN 
by IBonzini et al.l ()201Sf ). i.e. it must show radio emis¬ 
sion in excess of that produced by starbursts. This 
is measured by using the q 24 obs parameter, defined as 
g 24 obs = log(S' 24 obs/S'r), where S' 24 obs is the observed flux 
density at 24um and S r is th e observed flux density at 
1.4GHz. IBonzini et al.l (|2013f) define an object as radio 
loud if g 24 obs is 2cr off of the starburst locus, defined using 
the M82 template (see their Fig. 2). 

iv) We also impose the additional constraint that its 
X-ray luminosity must not exceed T 2 -iofeey = 10"^'* erg 
s“^. This is to ensure that the objects are consistent with 
typical low power radio galaxies at both low and high red - 
shifts fe.g. iBalmaverde et al.ll2005 iTundo et al.l 1201211 . 
The upper limit to the X -ray luminosity is al s o con - 
sistent with the results of iTerashima fc WilsonI (1200,81 1. 
These authors showed that the radio loud - radio quiet 
dichotomy can be redefined using the X-ray emission in¬ 
stead of the flux in the optical band. This is particularly 
useful for Type 2 AGNs such as the ones considered in 
this work, since the optical emission form the accretion 
disk is heayily obscured in these object. According to 
such a scheme, powerful radio loud AGNs are present 
for Rx > —3.5, where Rx is defined as the logarithm of 
the ratio between the radio luminosity {vh^) at 1.4GHz 
and the X-ray luminosity L 2 _io. At low radio powers, 
it is easy to confuse a powerful radio quiet AGN that 
shows some radio emission with a radio loud object that 
is intrinsically weak at all wayelengths. Note that our se¬ 
lection criteria for the radio and X-ray luminosities cor¬ 
respond to Rx > —2.3, which is a rather conseryatiye 
yalue. We prefer to follow a conseryatiye approach be¬ 
cause for low luminosity AGNs the transition between 
radio-quiet and radio-loud most likely occurs at higher 
yalues than for high power objects, similarly to the classi¬ 
cal radio-loudness paramete r R deriyed using the radio- 
to-optical luminosity ratio (ITerashima fc WilsonI 120031 : 
iChiaberge et al.l 120051 : iSikora et al.ll2007lb 

v) The source must be obseryed with HST WFG3-IR. 
Mos t of the images are taken from the CAND FLS sur- 
yey (iKoekemoer et al.ll2011l : IGrogin et al.ll20lill , or from 
surrounding fields obseryed as part of program GO-12866 
(see Table |2l). The ECDF sample includes 6 HzLLRGs. 

We also use the same criteria to select HzLLRG in the 
CANDELS UDS field in the Subaru/XMM Newton Deep 


field surv ey fSXDF. lUeda et al.1 1200811 . The radio data 
are from (iSimpson et al.ll2012D . the 24/4m Spitzer IR data 
are taken from the Spitzer Public Legacy Survey of the 
UKIDSS Ultra Deep Survey (SpUDS, P.I. J.S. Dunlop) 
using the IRSA database. Only two galaxies that satisfy 
our selection criteria lie in the region covered by the HST 
WFG3 observations. 

The properties of the 8 galaxies included in the 
HzLLRG sample are reported in Table [2j In Fig. [T] 
we show the radio power vs. redshift distribution of 
these AGNs, compared with the other radio loud sam¬ 
ples. Note that the two radio loud samples at z > 1 
used in this paper (i.e. the 3GR and the HzLLRGs) are 
separated by about 4dex in radio power, on average. 

2.3. Radio quiet AGNs: low-power Type 2 AGNs at 

1 < 2 < 2.5 

The sample of low-power Type 2 AGN (LPTy2AGN) 
is selected from the 4Msec GDFS catalog and source clas¬ 
sification (jXue et al.ll2011h . We include objects that are 
classified as AGNs, that do not show strong broad emis¬ 
sion lines in the optical spectrum, and whose redshift is 
1.0 < z < 2.5. We also use the additional constraint 
that the intrinsic (de-absorbed) X-ray luminosity (inte¬ 
grated between 2 and 10 keV, in the rest frame of the 
source) must be L 2 -iokeV < 2 x 10“^^ erg s“^, in order 
to be consistent with the properties of the corresponding 
radio loud sample described in Sect. 12.21 The L 2 _io/cey 
luminosity is derived from the 0.5-8keV luminosity listed 
in the catalog, converted to the 2-lOkeV band using a 
photon index P = 1.8. Furthermore, for the objects 
that are detected in the radio, we require that the X- 
ray radio loudness parameter Rx defined in Sect. 12.21 
is < —3.5. The sample includes 26 objects. These are 
very low power AGNs that are similar to Seyfert 2 galax¬ 
ies in the local Universe. The properties of these AGNs 
are reported in Tables [3] The HST/WFG3-IR images 
for this sample are taken from the GAND ELS survey 
(IKoekemoer et al.ll2011l : IGrogin et al.ll20lill . 

2.4. Radio quiet AGNs: high-power Type 2 AGNs at 

1 < z < 2.5 

As for the low power AGNs described above, the sam¬ 
ple of high-power Type 2 AGN (HPTy2AGN) is drawn 
from the GDFS. When selecting this sample we want to 
match as close as possible both the redshift range and the 
bolometric luminosity of the powerful 3GR radio galaxy 
sample. 

The sample of high-power Type 2 AGNs includes ob¬ 
jects classified as AGN with redshift 1 < z < 2.5 and 
intrinsic X-ray luminosity L 2 _iofcey > 10'^'^ erg s“^. 
The lower limit in X-ray power is chosen to match the 
properties of the 3GR sample described in Sect. 12.11 
Two of the objects (166 and 577) show relatively broad 
Mgll and GUI] emi ssion lines in the optical spectrum 
(|Szokolv et al.l[2003) . However, the HST images clearly 
show the host galaxy and there is no evidence for the 
presence of any strong unresolved nuclear source, as it 
would be expected for a powerful Type 1 QSO. There¬ 
fore, for the purpose of this work, we consider those two 
objects as Type 2, regardless of the presence of broad 
lines. 

As for the LPTy2AGN, for the objects that are de¬ 
tected in the radio, we require that the X-ray radio loud- 
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TABLE 1 

High z SCR, Radio Galaxies 


Name 

R.A. (2000.0) 
hh:mm:ss.ss 

Decl. (2000.0) 
dd:mm:ss.ss 

Redshift 

z 

log Pi. 4GHz 
log[erg Hz~^] 

^ 24fj,m 

log[erg Hz~^] 

3C 210 

08:58:10.0 

+27:50:52 

1.169 

35.13 

45.65 

3C 255 

11:19:25.2 

-03:02:52 

1.355 

35.29 

<44.54 

3C 257 

11:23:09.2 

+05:30:19 

2.474 

35.94 

45.94 

3C 297 

14:17:24.0 

-04:00:48 

1.406 

35.33 

44.84 

3G 300.1 

14:28:31.3 

-01:24:08 

1.159 

35.37 

- 

3G 305.1 

14:47:09.5 

+76:56:22 

1.132 

35.10 

45.36 

3C 322 

15:35:01.2 

+55:36:53 

1.168 

35.20 

44.90 

3C 324 

15:49:48.9 

+21:25:38 

1.206 

35.34 

45.48 

3G 326.1 

15:56:10.1 

+20:04:20 

1.825 

35.75 

45.64 

3C 356 

17:24:19.0 

+50:57:40 

1.079 

34.96“ 

45.52 

3G 454.1 

22:50:32.9 

+71:29:19 

1.841 

35.60 

45.67 


Note. — Si, 4 from IClondon et al.l lll99SI^ except for ° ILaing fc Peacocld J198(1L The 24/rm 
luminosities are derived from the fluxes published in IPodigachoskr^r'aLn'20151 1 . 


TABLE 2 

High z Low Luminosity Radio Galaxies 


ID 


R.A. (2000.0) Decl. (2000.0) Redshift HST Survey or Prog. ID log Pi. 4 GHz log L 2—10 q24 

hh:mm:ss.ss dd:mm:ss.ss z log[erg s“^ Hz“^] log[erg s"*^] 


EGDFS 


65 

03:31:23.30 

-27:49:05.80 

1.10“ 

12866 

31.85 

<42.56 

-1.54 

127 

03:31:34.13 

-27:55:44.40 

1.06 

12866 

30.51 

<42.69 

-0.09 

215 

03:31:50.74 

-27:53:52.15 

1.77 

12866 

31.13 

<42.69 

-0.46 

338 

03:32:10.79 

-27:46:27.80 

l.Ol*" 

GANDELS 

31.29 

42.43 

-0.98 

410 

03:32:19.30 

-27:52:19.38 

I.IO*" 

GANDELS 

30.37 

<42.27 

-0.23 

412 

03:32:19.51 

-27:52:17.69 

1.06 

GANDELS 

30.93 

<42.24 

-0.51 

UDS 

48 

02 18 18.38 

-05 15 45.2 

1.56 

CANDELS 

32.27 

<43.93 

0.09 

124 

02 17 04.77 

-05 15 18.1 

1.28*’ 

CANDELS 

31.54 

<43.73 

-0.50 


Note. — For the EGD FS sources the ID cor responds to the source ID in the lBonzini et aO II2012I1 ca talog. For the 
UPS g alaxies, the ID is the lSimoson et ^ II2012I 1 source number. Redshifts are photometric redshifts from IBonzini et al.l 
II2OI2IT . except where stated otherwise. The UDS sources are undetected in the X-rays. 

® Photometric redshift from ICardamone et a Denia) 

Spectroscopic redshift 


ness parameter Rx defined in Sect. 12.21 is < —3.5. This 
sample includes 9 AGNs (see TableS]). The HST/WFC3- 
IR ima ges for this sample are taken from the CAN DELS 
survey ( Koekemoer et al.ll2Qlll : iGrogin et al.ll201lD . 

Since the X-ray luminosity might be poorly con¬ 
strained because of the uncertainty on the amount of 
nuclear obscuration, we should also make sure that this 
sample and the Hz3Cs have similar mid-IR luminosities. 
We performed a K-S test to check whether the distribu¬ 
tions of 24/rm luminosity of the Hz3C and HPTy2AGN 
are different The resulting p-value is p=0.1, therefore 
we cannot reject the null hypothesis that the two distri¬ 
butions are drawn from the same population. 

2.5. Non-active galaxies at 1 < z < 2.5 

This sample is derived from the 3D-HST s urvey 
of the GOODS-SOUTH held (iGiayalisco et al.l 12004 
iBrammer et al.l l2012t iSkelton et al.l 1201^ We select 
galaxies with spectroscopic redshift 1 < z < 2.5 and 

^ Note that we did not perform the same test for the LPTy2AGN 
because those are all undetected at 24/im. 


with magnitude in the F140W filter 19 < mpiAow < 22. 
The magnitude range is chosen to match the magnitudes 
of the 3CR radio galaxies from Sect. 12.11 (Hilbert et ah, 
in prep.). The full sample includes 145 galaxies. We 
limit the sample to the 50 objects included in the south¬ 
ern half area of the held observed with WFC3-IR and 
the F140W hlter as part of the 3D-HST survey. Five 
of these objects are AGNs, therefore we are left with a 
sample of 45 galaxies. These are sufficient to provide 
us a statistically sound sample to be compared with the 
active galaxy samples described above. We refer to this 
sample as the bright sample of non-active galaxies. Note 
that we could in principle derive a larger sample by us¬ 
ing the full area covered by the SD-HST image of the 
GOODS-S held. However, this would not signihcantly 
improve the statistics for a sample that is already the 
largest we consider. 

In the same area of the sky, we also select a compara¬ 
ble sample of fainter galaxies (22 < itik^ < 24) from the 
same HST survey. This sample spans the same magni¬ 
tude range as the bulk of our radio-quiet AGN popula¬ 
tion. In the following, we will refer to this sample as the 
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TABLE 3 

Radio Quiet Low Power Ty2 AGN at 1 < x < 2.5 


ID R.A. Decl. Redshift log L 2 —lo 

hh:mm:ss.s dd:mm:ss.s z log[erg s~^] 


Low Power Type 2 AGN (GDFS) 


125 

03:32:06.77 

-27:49:14.10 

1.050 

42.04 

147 

03:32:09.22 

-27:51:43.50 

1.352 

41.99 

216 

03:32:15.26 

-27:44:38.60 

1.109 

41.70 

225 

03:32:15.91 

-27:48:02.20 

1.520 

41.90 

226 

03:32:16.04 

-27:48:59.90 

1.413 

41.91 

247 

03:32:17.84 

-27:52:10.80 

1.760 

42.15 

317 

03:32:23.16 

-27:45:55.00 

1.224 

42.25 

318 

03:32:23.17 

-27:44:41.60 

1.571 

42.09 

321 

03:32:23.61 

-27:46:01.40 

1.033 

41.98 

376 

03:32:27.04 

-27:53:18.60 

1.103 

41.88 

389 

03:32:28.62 

-27:45:57.20 

1.626 

41.97 

394 

03:32:28.85 

-27:47:56.00 

1.383 

41.89 

416 

03:32:29.94 

-27:52:52.80 

1.017 

41.76 

419 

03:32:30.05 

-27:50:26.80 

1.005 

41.48 

428 

03:32:31.11 

-27:49:40.00 

1.508 

41.79 

455 

03:32:33.06 

-27:48:07.80 

1.188 

41.75 

462 

03:32:33.67 

-27:47:51.20 

1.388 

41.84 

463 

03:32:33.84 

-27:46:00.60 

1.903 

42.09 

491 

03:32:35.80 

-27:47:35.10 

1.223 

42.13 

493 

03:32:35.98 

-27:48:50.70 

1.309 

42.19 

504 

03:32:36.35 

-27:49:33.40 

1.508 

41.81 

536 

03:32:39.07 

-27:53:14.80 

1.380 

41.98 

541 

03:32:39.42 

-27:53:12.70 

1.381 

42.00 

545 

03:32:39.65 

-27:47:09.60 

1.317 

41.88 

558 

03:32:41.01 

-27:51:53.40 

1.476 

42.19 

579 

03:32:43.45 

-27:49:02.20 

1.603 

41.97 


Note. — The ID corresponds to the source ID in the 
IXue et ^ II 2 OIII ') catalog. The redshifts (z^dopt in the 
Xue et al. catalog) are spectroscopic, if available, or pho¬ 
tometric. The intrinsic X-ray luminosities (converted to 
the 2-lOkeV ban d as explained in the text) are also from 
IXue et ^ 1 I 2 OIIII . 


faint galaxies sample. 

Even if we are selecting objects in the very same red- 
shift range, the near-IR (rest-frame optical) magnitude 
might not be a good tracer of the stellar mass because 
of the possible presence of obscuration in some objects. 
We do not believe that this might significantly affect our 
results. However, in order to perform a sanity check, 
we also select a sample of non-active galaxies in the 
full GOODS-S area of 3D-HST matched to the stellar 
mass estimates of the Hz3C galaxies. At 1 < z < 2.5 
these are typically b etween 1 x 10^^ and 5 x 10^^ Mq 
(|S evmour et al.l[200^ . although a small number of lower 
mass objects are also present. Although stellar mass 
measurements based on SED fitting heavily rely on mod¬ 
els, it is important to check that our non-active galaxy 
samples are correctly representing the population of ob¬ 
jects we need for comparison. If we restrict our range 
of masses between the above values, and after excluding 
those that are AGNs, we obtain 23 objects. About half of 
them are in common with the bright sample. While we do 
not consider this as our main control sample of non-active 
galaxies, results for this sample are briefly discussed in 
Sect, m for the sake of completeness. 

The HST/WFG3-IR images are taken as part of the 
3DHST survey with the F140W filter. 

3. MERGER FRAGTION 


In order to measure the merger fraction we use at least 
four human classifiers for each sample. For the Hz3G 
and the Ty2 AGN samples we used six classifiers, to be 
sure that the authors of the paper were not biased a 
priori towards any specific result. Each classifier visually 
inspects all of the targets. While being qualitative in 
nature, visual classification has been proven to be an 
effective way of classify mergers, since the eye can pick 
all of the different signatures of such events. On the other 
hand, quantitative methods using the Gini coefficient G, 
the concentration index C, the asymmetry index A, or 
the second-order moment of the brightest 20 per cent of 
the light (M 20 ) only select particular merger s that each 
of th e above indicators are able to identify (jLotz et al.l 
[Ml . Therefore, if we use any of those methods we 
would derive smaller merger fractions. 

We classify objects as mergers (or post-mergers) if clear 
signatures of a mergers are present. These include the 
presence of double/multiple nuclei, close pairs, tidal tails, 
bridges, or distorted morphologies clearly indicating a 
recent or ongoing merger. Close pairs are defined using 
a projected separation of less than 25kpc between the 
center of the galaxies involved, corresponding to 3 arcsec 
in the redshift range spanned by our z > 1 samplefl 
Such a scale is similar to the typical s eparations observed 
in galaxy pairs at low redshifts (e.g. iSmith fc Heckm^ 
119891 iBehroozi et al.ll2015lPI . 

When possible, mergers are also visually distinguished 
between major and minor, assuming the usual separa- 
tion at a 4: 1 mass/apparent brightness ratio (e.g. 
iLotz et nil2011|) . If there is no clear evidence for any 
of the above properties, we classify the object as a non¬ 
merger. For the AGN samples, a blind classification is 
performed, i.e. classifiers did not know whether each ob¬ 
ject was radio-loud or radio-quiet, and if it belonged to 
the high or low-power class. Classifiers are also asked to 
classify objects more than once (typically twice). 

While substantial agreement exists between the differ¬ 
ent classifiers for most of the objects, different people 
may see different features in each image, thus the classi¬ 
fication for single sources may differ. This is why we do 
not report the merger classification for each galaxy in the 
tables. In order to reduce the effects of large deviations 
among the different classifiers we calculate 10%-trimmed 
means for each samnlJ^. We report the results for each 
of the samples at z > 1 in Tab. O 

® Between z=l and z=2.5 the projected scale varies by less than 
0.1 arcsec. 

® Note that the definitions of “double nucleus” and “close pair” 
are only formally different. For the purposes of our work, these 
objects are in fact the same. The different names only identify 
the appearance of the object in the image, since our classification 
is based on visual inspection. We classify the object as a double 
(or multiple) nucleus if the two objects are not clearly separated 
and they appear to lie within a common envelope. If the galaxies 
are well separated, but their nuclei are less than ~ 25kpc apart, 
then we call it a close pair. However, it is important to note that 
all of these objects are required to display evidence of bridges or 
asymmetries to be classified as mergers. 

Trimmed (or truncated) means are robust estimators of cen¬ 
tral tendency, and it is less sensitive on the outliers than the mean. 
Trimmed means are derived by calculating the mean after discard¬ 
ing parts at the high and low ends of a probability distribution. In 
our case, the distribution is defined by the number of mergers ob¬ 
tained for each sample by each of the classifiers, and we reject 10% 
at both ends. Note that a 50% trimmed mean would correspond 
to the median. 
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TABLE 4 

Radio Quiet High Power Ty2 AGN at 1 < z < 2.5 


ID 

R.A. 

hh.mm.ss.s 

Decl. 

dd.mm.ss.s 

Redshift 

z 

log L 2_10 
log[erg s-l] 

log L24^tm 
log[erg s-l] 

166 

03:32:10.93 

-27:44:15.20 

1.605 

44.21 

44.81 

243 

03:32:17.19 

-27:52:21.00 

1.097 

44.05 

44.39 

257 

03:32:18.35 

-27:50:55.61 

1.536 

44.07 

44.30 

278 

03:32:20.07 

-27:44:47.51 

1.897 

44.06 

45.04 

351 

03:32:25.70 

-27:43:06.00 

2.291 

44.31 

45.46 

518 

03:32:37.77 

-27:52:12.61 

1.603 

44.23 

45.38 

577 

03:32:43.24 

-27:49:14.50 

1.920 

44.12 

44.70 

713 

03:33:05.90 

-27:46:50.70 

2.202 

44.02 

- 

720 

03:33:07.64 

-27:51:27.30 

1.609 

44.39 



Note. — The ID corresponds to the source ID in the IXue et al.l 112011]') 
catalog. The redshifts (zadopt in the Xue et al. catalog) are spectroscopic, 
if available, or photometric. The intrinsic X-ray luminosities (converted 
to the 2-lOkeV band as explained in the text) are also from jXue et al.| 
goU). The 24/im luminosities are from fUardamone et al.l (l20Q8ll . 



Fig. 2.— RGB images of four radio galaxies from the high-z SCR sample showing clear evidence for recent or ongoing major merger. The 
objects are (from left to right, top to bottom) 3C 210, 3C297, 3C356 and 3C454.1. The HST WFC3-IR F140W image was used for the R 
channel. The WFC3-UVIS F606W image was used for both the G and B channels. North is up, East is left. Data from Hilbert et al. (in 
preparation) 


Examples of the morphologies observed in the differ- 
ent samples are given in Figs. [2l [3l and 21 In order to 
avoid confusion between our merger classification crite¬ 
ria, we use the figures to give a few specific examples. 
3C 297 and 3C356 in Fig. [21 are both classified as merg¬ 
ers based on those clear signatures of interaction. 3C 297 
also shows a double nucleus. Although not all cases are 
as straightforward as these, similar considerations can be 
made for objects in the other samples (see e.g. the sys¬ 
tems shown in the right panes of both Figs. [3] and HI)- In 
particular, the objects shown in Fig. [3] are most likely ex¬ 
amples of dry mergers in which the galaxy isophotes are 
highly asymmetric (see e.g. lBell et al.ll200m . The central 
panels of Fig. 21 show objects that some of the classifiers 
classified as mergers, while others did not. Those that 


classified the objects as mergers noticed some asymme¬ 
tries in the isophotes of the galaxies, in addition to the 
presence of small companions. 

The observed merger fraction in the radio-loud samples 
(Hz3C and HzLLRG) are clearly larger than those found 
for radio-quiet and non-active galaxies. In Sect. 21 'we 
test this result using careful statistical analysis. 

3.1. On the impact of different sensitivities on the 
merger classification 

The goal of this work is to identify merging systems 
in samples of objects that were observed as part of dif¬ 
ferent surveys or observing programs. All of the images 
were taken using the WFC3-IR camera. Its extremely 
high sensitivity and spatial resolution allows us to detect 
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Fig. 3.— Four examples of the morphologies observed among 
the 2 > 1 radio quiet AGNs. In the top panels, we show two 
LPTy2AGNs. A non-merger and a merger are shown on the left 
and on the right, respectively. In the bottom panels, the same is 
shown for two HPTy2AGNs. The circles are 6 arcsec radius, which 
correspond to ~ 50 kpc at the redshift of the objects. Images 
are from CANDELS IIKoekemoer et ^120111 ; IGrogin et al.ir2011l . 
HST/WFG3/F160W). 


TABLE 5 

Observed Merger Fractions (trimmed means) for 
THE 2 > 1 samples 


Sample 

Sample Size 

Merger Fraction 

Hz3G 

11 

100% 

HzLLRG 

8 

88% 

LPTy2AGN 

26 

38% 

HPTy2AGN 

9 

33% 

Bright Galaxies 

45 

27% 

Faint Galaxies 

46 

20% 

RL (3G -1- LLRG) 

19 

95% 

RQ (HP-|-LPTy2AGN) 

35 

37% 


low surface brightness features that characterize recent 
merger events at z ~ 1 — 2 even in short (1 orbit or less) 
observations. However, for the sake of clarity, we report 
in Table |6] the 5cr surface brightness limits for the dif¬ 
ferent surveys we use. To derive the limits we used the 
WFC3-IR Exposure Time Calculator (ETC). We assume 
a 2x2 pixel extraction area and a spectrum of an ellipti¬ 
cal galaxy with a UV upturn to perform the conversion 
between E160W (used for the CANDELS observations) 
and F140W magnitudes. Although the CANDELS deep 
images we used for the radio quiet samples and for the 
HzLLRGs are deeper than all other data, we checked 
that the merger classification is not different if the shal¬ 
lower 3D-HST images are considered. The images show 
the very same features, irrespective of the redshift of the 
source. The faintest tidal structures are less prominent 
in the shallower images, but that does not significantly 
affect our classification. The short exposure times of the 
3C snapshot survey are also sufficient to detect the faint 
features we are interested in, in the range of redshift of 
our sample, as shown by the fact that that sample has 
the greater observed merger fraction. 

4. STATISTIGAL ANALYSIS 

The main goal of this work is to investigate if mergers 
are associated with AGN activity and if that plays a role 
in triggering such a phenomenon. Most importantly, we 


want to test whether different classes of AGNs (e.g. radio 
quiet vs. radio loud, low vs high bolometric luminosity) 
are more likely to be triggered by merger than others. 
We also test the hypothesis that AGNs are no different 
than non-active galaxies. 

We perform a set of statistical tests to compare the 
derived merger fractions for the five different samples 
described in Sect. [2] Throughout the paper, statistical 
tests are perf ormed using differe nt techniques in the R 
environment (jB Core TeamllMT^ . 

We use the Bayesian version of the proportion test us¬ 
ing bayes.prop. test, as par t of the Bayesian First Aid 
package for R (|Baathll2014D . In principle, Bayesian tests 
are more useful than classical proportion tests, since they 
provide a n estimate of the relative frequency of suc¬ 
cess (e.g. iBolstadI I2007D . In this simplified version of 
the Bayesian tests, the priors are uninformative, i.e. a 
uniform distribution. This is suitable for our purposes, 
since we have no a priori knowledge of the distribution 
of mergers in each samplJ*~*1. 

In table [7] we list the relative frequencies of success for 
each sample, together with the 95% credible intervals. 

By testing each sample against each of the others, we 
can firmly establish that the Hz3C has a larger merger 
fraction than both of the Ty2AGN samples (P> 99.9%). 
The Hz3C merger fraction is also larger than that of 
the non-active galaxies (both bright and faint samples, 
P> 99.9%). Eurthermore, the HzLLRGs are significantly 
more associated with mergers than the non-active galax¬ 
ies (P> 99.9). All other tests are inconclusiv43. How¬ 
ever, these are very important results. 

Eurthermore, if we merge the two samples of z > 1 
radio loud objects and we test them against the radio 
quiet AGNs, the result is extremely robust. The merger 
fraction in our sample of radio-loud AGNs is significantly 
higher than that in the radio-quiet sample. This is again 
a notable result, since the (Bayesian) merger fractions 
for the RL and RQ sample are 92% and 38%, respec¬ 
tively. As expected, the same result holds (with an even 
higher statistical significance) for the RL sample against 
the non-active (both bright and faint) galaxies, while 
the RQ sample does not show any statistically significant 
difference with respect to the non-active galaxies. It is 
particularly important to note the results for the samples 
that are matched in magnitude. On the one hand, the 
RL and the bright galaxy sample are different, and on 
the other hand, the RQ and the faint galaxy sample are 
statistically indistinguishable. 

As a sanity check, to avoid possible biases from a se¬ 
lection made using near-IR magnitudes, we also used 
the smaller sample of high-mass non-active galaxies (see 
Sect. 12.5]) . In particular, given that we measure a merger 

We also perform a complete set of classical proportion tests 
using the R task prop.test. For these classical tests we reject the 
null hypothesis that the merger fractions in two samples are the 
same if the p-value is p < 0.01. Not surprisingly, the results are 
perfectly in agreement with the Bayesian relative frequencies. For 
the sake of clarity, here we only discuss in details the Bayesian 
results. 

For a 2-sample proportion test in classical statistics this means 
that we cannot reject the null hypothesis that the observed merger 
fractions in two samples are the same (p-value < 0.01). For our 
Bayesian treatment, this implies that the probability that one sam¬ 
ple has a higher (or lower) merger fraction than the other is less 
than 99%. 
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Fig. 4.— Six examples of non-active galaxies from our z > 1 sample. Both in the top and bottom rows, objects that are unanimously 
classified as non-mergers or mergers are shown in the left and right panels, respectively. Objects with mixed classification are shown in 
the middle panels. The circles are 6 arcsec radius, which corresp ond to ~ 50 kpc at the redshift of the objects. Images are from 3D-HST 
IIGiavalisco et ani2004l : iBrammer et al.ll2012 : [Skelton et al.ll2od HST/WFC3/F140W) 


TABLE 6 

Surface brightness limits 


Survey or HST Prog. ID 

Sample 

HST Camera/Filter 

5(T Surface Brightness Limit 
MF140W [ABmag arcsec"^] 

Sensitivity 
AB mag 

SNAP-13023 

Hz3C 

WFC3/F140W 

23.9 

25.8 

CANDELS wide 

KQ and HzLLHG 

WFC3/F160W 

24.1 

26.5 

CANDELS deep 

RQ and HzLLRG 

WFC3/F160W 

25.6 

27.2 

3D-HST 

non-active 

WFC3/F140W 

24.1 

26.2 

GO-12866 

HzLLRG 

WFC3/F160W 

25.0 

26.7 

Low-z Samples 

GO-9045 

3CRR,6C,7C,TOOT 

WFPC2/F785LP 

20.7 

24.9 

SNAP-10173 

SCR (z< 0.3) 

NICMOS/NIC2/F160W 

22.5 

24.3 


Note. — The 5cr limits are estimated using the WFC3-IR ETC and assuming a 2x2 pixel extraction region. In column 
4 we report the limit surface brightness for each datasets, converted to the WFC3-IR F140W filter to allow for an easy 
comparison between the different surveys. An elliptical galaxy spectrum with UV upturn redshifted to the appropriate 
redshift for each sample was used for the conversion between the two WFC3 filters. RQ in the samples column refers to 
both the LP and HPTy2AGN samples. For the low redshift samples (GO9045 and SNAP10173) we used synphot to convert 
magnitudes to the WFC3 filter system. In column 5 we report the 5 (t image sensitivities calculated for point sources in 
each of the bands used for the observations. 


fraction of the high-mass sample is statisti¬ 

cally different from the Hz3C, with a Bayesian probabil¬ 
ity > 99.9%. On the other hand, the merger fraction 
in such a sample is statistically indistinguishable from 
those in both the faint and the bright samples of non¬ 
active galaxies. 

Finally, we wish to point out that among the sam¬ 
ples of radio-loud AGNs at z > 1 the large majority of 
the objects 90%) appear to be associated with major 
mergers. On the other hand, for radio quiet AGNs at 
z > 1, only ^ 50% of the observed mergers are major 
mergers. Note that this is only based on our qualitative 
visual classification of the mergers. A more careful clas¬ 


sification would require 2-D galaxy modeling and SEDs 
to derive stellar mass estimates for the galaxies involved 
in each merger, which is beyond the goals of this work. 

4.1. Lower redshift samples of radio loud AGNs 

One important goal of this work is to establish whether 
the merger fraction for radio-loud AGN depends on ei¬ 
ther redshift or luminosity. In order to test if that is 
the case, we used additional lower redshift samples, for 
c omparison with the z > 1 data. 

iMcLure et al.l (|2004f) observed a sample of 41 inter¬ 
mediate (z~ 0.5) radio galaxies with HST and WFPG2 
(GO 9045, P.I. Willott). The sample is taken from four 
different complete catalogs, spanning about 4 dex in ra- 
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TABLE 7 

Estimated Bayesian probabilities eor the z > 1 samples. 


Sample 

Merger Fraction 

95% Credible Interval 

Hz3C 

94% 

0.78 - 1.00 

HzLLRG 

82% 

0.57 - 0.99 

LPTy2AGN 

39% 

0.22 - 0.57 

HPTy2AGN 

36% 

0.11 - 0.64 

Bright Galaxies 

27% 

0.15 - 0.40 

Faint Galaxies 

20% 

0.10 - 0.33 

RL (3C -1- LLRG) 

92% 

0.78 - 1.00 

RQ (HP-|-LPTy2AGN) 

38% 

0.23 - 0.54 


TABLE 8 

Estimated Bayesian probabilities eor the z < 1 radio 
LOUD AGN samples. 


Sample Size Merger Fraction 95% Credible Interval 


Willott Samples 

3GRR 

13 

74% 

0.53 - 0.94 

6GE 

7 

68% 

0.38 - 0.93 

7GRS 

9 

65% 

0.37 - 0.89 

TOOT 

12 

65% 

0.39 - 0.87 

Full Sample 

41 

70% 

0.56 - 0.83 

Low Redshift Radio galaxies 

2 < 0.3 3GR 

101 

88% 

0.81 - 0.94 


dio power. The four catalogs are all complete, flux lim¬ 
ited and radio selected at low frequencies. The sample 
includes objects from the 3CRR, 6CE, 7CRS and the 
TexOx-1000 (TOOT sample, hereinafter) (|McLure et al.l 
12004 and references therein). The redshift range 
spanned by the objects observed with WFPC2 is between 
z = 0.40 and 0.59. We will refer to these galaxies as the 
Willott sample throughout the paper. 

For these objects, we retrieved the HST data from 
the MAST and we reduced them using Astrodrizzle 
(|Fruchter et al.ll2M^ . We then classify the objects be¬ 
tween mergers and non-mergers, since n o classification 
was p rovided in the original paper by iMcLure et al.l 
(|2004ll . and for consistency with the other samples. Note 
that the filter used for these observations (F785LP) pro¬ 
vides images at ~ 5300A in the rest frame. This is 
similar to the rest frame wavelength of the WFC3-IR 
observations for the z > 1 samples described above. In 
Tab. [5] we report the merger fractions for this sample. 
The merger fractions are statistically compatible with 
those for the radio-loud samples at higher redshifts. 

At even lower redshifts, we use the observations 
of the SCR sample with z < 0.3 taken with HST- 
NICMOS (prograrn SNAP10173) and the F160W filter 
(jMadrid et al.ll200'^ . Although these observations are in 
the rest-frame IR, they are better suitable for our pur¬ 
poses tha n the optical data of the same sample ta ken with 
WFPC2 (Ide Koff et al.lll99^IMartel et al.lll999[) because 
the NICMOS images are significantly deep er. For this 
sampl e, we consider the results published bv iFlovd et al.l 
(|2008f) . These authors classified all of the objects in 
mergers, pre-mergers, tidal-tails, major and minor com¬ 
panions. The field of view covered by NIC2 is 19.2" x 


19.2", which corresponds to 50x50 kpc^ at the median 
redshift of the sample. They found that 89 out of 101 ob¬ 
jects fall into at least one of these categories. This corre¬ 
sponds to a merger fraction of 0.88, with a Bayesian 95% 
credible interval [0.81-0.94]. If we exclude objects with 
only minor companions (i.e. > 1 mag fainter than the 
radio galaxy host, as defined by Floyd at ah, 2008), the 
fraction is reduced to 0.76 with a Bayesian 95% credible 
interval [0.68-0.84]. 

Since the observations are not homogeneous, these two 
lower redshift samples are only considered here for com¬ 
parison. However, we do not base our conclusions on 
these samples only. 


5. DISCUSSION 

The main result of this work is that the samples of ra¬ 
dio quiet and radio loud Type 2 AGNs at 1 < z < 2.5 
have different merger fractions. We showed that there is 
clear statistical evidence that the radio loud AGNs al¬ 
most always reside in environments where mergers are 
undergoing, or that recently happened. As discussed 
below, this is independent of the radio (or bolometric) 
power of the AGN, and it is statistically compatible with 
merger fractions as high as ^ 100%. 

The sample with the highest observed merger fraction 
is the z > 1 3GR. In principle, this sample could be bi¬ 
ased, since HST only imaged ~ 35% of the radio galaxies 
with z > 1 included in the 3GR catalog. However, it is 
important to note that the observations were performed 
as part of an HST SNAPSHOT program. In SNAP pro¬ 
grams, targets are randomly selected based on the avail¬ 
ability of gaps in the HST schedule. Therefore, there 
was no a priori knowledge of the properties of the ob¬ 
served targets with respect to the complete sample. But 
since the observed sample is small, it is still possible that 
we ended up picking a large number of objects in merg¬ 
ers only by chance. In the following we briefly discuss 
such an issue. For example, we could in principle as¬ 
sume that the original population is composed by ^50% 
of mergers and 50% of non-mergers, similarly to what is 
observed in our sample of non-active galaxies. We can 
test the probability of obtaining a sample of 12 merg¬ 
ers, randomly extracted from the complete sample of 34 
radio galaxies included in the 3CR catalog with z > 1. 
The selection of the 3GR catalog only covered part of 
the sky (Deck > —10 deg). Therefore, we should cor¬ 
rect for the area coverage to obtain the total number 
of radio sources (i.e. 58 radio galaxies). The binomial 
probability of obtaining 12 mergers out of 12 observed is 
extremely small (P = 2 x 10“^). However, when objects 
are randomly drawn from a sample without replacement, 
the Hypergeometric distribution must be used instead of 
the binomial. In that case, the probability of obtain¬ 
ing 12 mergers out of 12 observed, from a population of 
58 objects in which 50% of the objects are mergers is 
P = 5.8 X 10“®. Note that if the merger fraction in the 
complete sample is higher than 72%, then the (hypergeo¬ 
metric) probability of obtaining 12 mergers is P > 0.01. 
Since that corresponds to the significance level we set 
for all of the statistical tests, we can state that we can¬ 
not completely rule out that the merger fraction in the 
3CR is as low as ~ 72%. Interestingly, this is similar 
to the lower value of the credible interval given by the 
Bayesian analysis for this sample (see Tab.]?]). The only 
































RLAGN are mergers 


11 


other selec t ion bi as that might affect our sample is the 
lEddingt^ (jl913ll bias. However, that would have the 
effect of lowering the actual number of sources in our 
sample, and that in turn would lower the probability of 
the hypergeometric distribution. 

Therefore, we conclude that even if the number of SCR 
galaxies observed at z > 1 is small, the probability that 
the observed merger fraction is overestimated because of 
any selection bias is extremely small. 

5.1. Radio loud samples: no trends with redshift and 
luminosity 

In this section we discuss the results we obtained for 
the different radio loud samples. It is in fact particu¬ 
larly interesting to investigate whether the merger frac¬ 
tion in these object may depend on e.g. redshift, lumi¬ 
nosity, or on the original criteria for each sample selec¬ 
tion. However, while the statistical analysis clearly shows 
that radio loud objects are almost all associated with 
mergers, one important caveat is that the samples are 
small. This is particularly true for the radio-loud sam¬ 
ples (except for the low-z SCR), but also for the high-z 
radio-quiet comparison sample. Unfortunately, this pre¬ 
vents us from drawing statistically firm conclusions on 
any possible trends between the different samples of ra¬ 
dio galaxies. For example, we cannot determine whether 
the HzLLRC’s are less likely to be associated with merg¬ 
ers than the HzSC’s even if the observed fraction is lower 
in the former than in the latter. Similarly, we cannot 
determine whether each of the samples at z ~ 0.5 that 
compose the Willott sample behave differently, for exam¬ 
ple as a result of their different radio power. 

However, it is interesting to consider radio-loud sam¬ 
ples grouped by redshift bin. We can thus compare the 
merger fractions for the low-z SCR, the whole Willott 
sample (TOOT-I-7CRS-I-6CE-I-3CRR), and the high-z 
objects (3CR-l-HzLLRC). In Sect. 0] we showed that we 
find no statistical evidence that that the observed merger 
fractions are different for any of those groups. Therefore, 
while we wish to stress that this does not imply that they 
are the same, we can conclude that the data show no ev¬ 
idence for a redshift evolution. By comparing the 95% 
Bayesian credible intervals, we can also conclude that 
the merger fractions do not differ by more than ~ 20%. 
Since the samples are well separated in radio luminos¬ 
ity, it is straightforward to perform a similar analysis for 
samples grouped by radio luminosity, and conclude that 
the merger fraction in low and high-power samples does 
not differ by more than 20%. This is a notable result, 
since these samples are separated by more than 4 dex 
in radio power, and span a wide range of redshift (see 
Figll]). 

In Figl5] we show the Bayesian 95% credible inter¬ 
vals for each of the groups of radio loud objects plot¬ 
ted against the redshift (left panel) and radio luminos¬ 
ity (right panel) range spanned by each group. In the 
left panel (merger fraction vs redshift) we also include 
the radio-quiet ACNs and the two samples of non-active 
galaxies. The radio-quiet ACNs and the non-active 
galaxies are not plotted in the right panel, since any 
trends with the radio power would be irrelevant. In fact, 
the origin of radio emission in radio-quiet ACN is still de¬ 
bated, as it could be either thermal or non-thermal, and 
a possible contribution from starbursts could not be ex¬ 


cluded at the lowest luminosities. Starbursts are instead 
the most likely origin for the radio emission in non-active 
galaxies. These objects would lie on the bottom-left of 
the figure in the right panel, but any correlation with 
the merger fraction would be meaningless. As it is clear 
from the figure (left panel), the only group of radio loud 
ACNs that is still marginally compatible with the ra¬ 
dio quiet samples is the Willott sample. However, this 
might be due to the fact that the images were taken 
with WFPC2, which was significantly less sensitive than 
WFC3-IR. In Tab. [5] we report the limit magnitude es¬ 
timated using the WFPC2 Exposure Time Calculator, 
and converted to the FI40W filter using synphot to al¬ 
low an easier comparison with the WFC3 observations. 
The value mFi 40 u; = 20.7 mag arcsec“^ was derived us¬ 
ing an elliptical galaxy spectrum redshifted to z=0.5. As 
a result of that, some images might not show faint sur¬ 
face brightness structures such as e.g. asymmetries or 
tidal tails, which might lead us to misclassify some of 
the objects as non-mergers. Deeper images with WFC3 
or with ACS in the I band should be taken in order to 
achieve a more reliable estimate of the merger fraction 
in the Willott sample. 

The right panel of Fig. [S] shows the merger fraction for 
each sample, against radio power at I5IMHJ*^. Here we 
plot the Willott samples separately, since they belong to 
different luminosity bins. The low-z 3CR span a large 
range in radio power. Most of these SCR objects are 
confined to the lower luminosity bin, as it is clear from 
Fig.m As noted above, except for the low-z SCR for 
which the merger fraction is much better constrained, 
the number statistics is small and the error on the merger 
fraction is large. However the figure clearly shows that 
all of the samples are located in the upper part up the 
diagram, and that no trend with radio power is visible. 

Finally, in Fig. [6] we show the merger fraction against 
the radio loudness parameter Rjf. In order to represent 
the average value of R^ for each sample, we calculated 
average values for both the radio and the X-ray lumi¬ 
nosities (see the caption for references). The data are 
very uncertain, since the information is extremely sparse. 
This especially holds for the X-ray data of the Hz3C and 
the Willott sample, while most of the HzLLRC only have 
upper limits in the X-ray band, so the radio loudness of 
the HzLLRCs is represented as a lower limit. Note that 
^ 20% of the radio-quiet objects are detected in the ra¬ 
dio band. The value of R^, for the radio-quiet samples is 
thus calculated using the average radio luminosity for the 
detected objects, and the points in the figure are shown 
as upper limits for the radio loudness parameter. For 
all other samples, to be safe, we assume uncertainties up 
to ~ ±0.5dex in R^. This is reasonable considering the 
uncertainties in the X-ray luminosity and the range of 
radio and X-ray luminosity spanned by the objects. 

This plot summarizes the main result of this work, 
i.e. that radio quiet sources are systematically associ¬ 
ated to smaller merger fractions (and they are located in 
the lower-left side of the figure), while radio loud ACN 
are unambiguously associated with mergers (and they are 
located in the upper-right side of the figure). 

Measurements at radio frequencies different form 151MHz 
were converted to the reference wavelength using a spectral index 
a = 0.8 and Fu oc 
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Fig. 5.— Left: merger fraction vs. redshift for the samples of radio loud and radio quiet AGN, and for non-active galaxies. The 
coordinates of the filled dots are the median redshift of each sample, and the Bayesian merger probability. The errorbars show the range of 
redshift spanned by each sample and the Bayesian 95% credible intervals for the merger probability. The High-z radio-loud sample (green 
dot) corresponds to the SCR at 2 : > 1 merged with the HzLLRG sample to improve the statistics. The radio-quiet Ty2 AGN samples at 
2 : > 1 are also merged (light red). The bright and faint samples of non active galaxies are plotted in yellow and cyan. The Willott sample 
(SCR, 6CE, 7CRS and TOOT samples at 2 : 0.5) is in red. The blue dot represents the SCR sample with 2 : < O.S. Right: merger fraction 

vs. radio luminosity at 151MHz for the radio loud samples. In this panel, the Willott samples (here marked as TOOT, 6C, 7C and SC) 
and the high-z radio-loud samples (HzSCR and HzLLRG) are plotted separately. The color code is the same as for the left panel. 
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Fig. 6.— Merger fraction vs. average radio loudness parameter 
Rx for the different AGN samples. Filled symbols are the radio 
loud samples, empty symbols are radio quiet. The Hz3C sam¬ 
ple is plotted as a circle, the HzLLRG is the square (lower limit), 
Willott’s 3GRR is the triangle, and the low-z 3CR is the pentagon. 
For the radio quiet samples (empty symbols) the LPTy2AGN 
sample is the triangle, and the HPTy2AGN is the square. The 
dashed line represents the radio loudness threshold for PG QSOs 
IITerashima fc Wilson|[2^003IL The solid line marks the 60% merger 
fraction that appears to roughly separate radio-loud and radio¬ 
quiet samples. Color code is the same as for Fig.[5] 


5.2. Comparison with other recent results 

Our results basica l ly ag ree with the findings of 
lR,amos Almeida et al.l (j2012ri and iRamos Almeida et al.l 
(|201,‘4) . Those studies focus on samples of relatively low- 
redshift {z < 0.7) sources with deep ground based ob¬ 
servations. These authors found that the large major¬ 
ity 80%) of radio-loud objects show disturbed mor¬ 
phologies, and they also reside in dense environments. 
However, the parameter space covered by those papers is 
limited to luminous low-redshift objects, while our work 
spans a significantly larger range both in redshift and 
luminosity. The s ame group also stud ied a sample of 
20 Type 2 quasars (iBessiere et al.ll201^ that includes at 
least two RLQS Os, based on th eir location in the Jj 5 GHz 
V L[o///] plane (|Xu et al.lll999D . As expected in light of 
our results, both RLQSOs are associated with mergers. 
The merger fraction in the full sample is as high as 75%. 
If we only exclude those two sources that exceed their 
definition of radio-loudness, the observed merger frac¬ 
tion among the remaining radio-quiet objects is still quite 
large (72%). While this might appear in disagreement 
with our measured fraction for the RQAGNs, the uncer¬ 
tainty on the value found by these authors is large, as a 
result of the small number statistics (13 possible mergers 
out of 18 objects). Since the samples are small, we can¬ 
not reject the null hypothesis that a 72% merger fraction 
for that sample is different from the fraction measured 
in our sample of RQAGNs. 

One of the major points here is that we showed that 
there is no evidence for a trend with luminosity, at least 
for the radio loud samples. This is ap parently at odds 
with the results of iTreister et al.l (I2012D . where only the 
most powerful samples appear to predominantly reside 
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in merging systems. However, it is interesting to briefly 
discuss their work in the light of our results. Two sam¬ 
ples show a significantly greater merger fraction than any 
other samples treated in that work. Those are the s ample 
of dust- reddened Type 1 QS Os of lUrrutia et al.l (120081 1 
and the iBahcall et al.l (jlOOTfl sample of Type 1 QSOs. 
The former is most likely biased in favor of a high merger 
fraction, because of the nature of those obscured quasars. 
The latter, instead, is more interesting. Of the 20 QSOs 
observed with HST bv IBahcall et al.l (jlOOTl l. 14 are radio 
quiet and 6 are radio loud. All of the 6 RLQSOs are ap¬ 
parently merging or show irregular features that might 
be explained with a merger, in agreement with our re¬ 
sults. Thus if we only limit the analysis to the RQQSOs, 
the fraction of mergers is ^ 50%, in agreement with our 
hndings. Therefore, we argue that the high fraction of 
mergers in the full Bahcall sample (65%) could in princi¬ 
ple be explained, in light of our results, by the fact that 
a significant fraction of the objects are radio loud. 

6. IMPLICATIONS FOR THE ORIGIN OF 
RADIO-LOUDNESS 

Our results show that radio galaxies are unambigu¬ 
ously associated with mergers, independently of redshift 
and luminosity, while radio-quiet AGN at z > 1 are in¬ 
distinguishable from normal galaxies in the same red- 
shift range. This result may have profound implications 
for our understanding of the mechanisms that enable the 
production of powerful relativistic jets from supermassive 
BHs. In fact, such a strong connection between mergers 
and radio loud AGNs may be a clue for a direct link 
between these two phenomena. 

The central question here is whether mergers may pro¬ 
vide, or at least substantially contribute to the physical 
conditions that ultimately enable the formation of jets in 
RLAGNs. 

6.1. Not all mergers may generate a RLAGN 

First of all, we should point out that most mergers do 
not generate a RLAGN. This is clear from the fact that 
a fraction of non-active galaxies at z > 1 are seen to be 
merging, but they show no signs of radio-loud activity. If 
there is an association between these two phenomena, it 
is not a univocal cause/effect relationship. Based on our 
results, the same also holds for radio-quiet AGNs since a 
fraction of those are associated with mergers. Thus, we 
conclude that mergers are unrelated to radio-quiet AGNs 
or, alternatively, only a fraction of those may be triggered 
by mergers. However, it is worth noting that our defini¬ 
tion of merger includes objects that still have to merge 
as well as objects for which the signs of a past merger 
are somehow still visible. Therefore, the timescale dur¬ 
ing which we observe a merger is probably of t he order of 
at least ~ 1-2 Gyr (e.g iDi Matteo et '^1200511 . The time 
scale for radio loud activity is most likely significantly 
shorter (~ 10^ — 10® yr) . Therefore, we cannot exclude 
that some of the non-active galaxies that we observe in a 
merger phase are turned-off radio loud AGNs, or, alter¬ 
natively, they still have to be turned-on. Summarizing, 
we believe that not all mergers may directly generate a 
radio-loud AGN. Below we discuss a few possible condi¬ 
tions for a merger to trigger RLAGN activity. 

6.2. Conditions for mergers to trigger RLAGNs 


We argue that when certain conditions are met, merg¬ 
ers may trigger radio loud nuclear activity. What we ulti¬ 
mately wish to know is what those conditions are. An im¬ 
portant piece of information here is that the association 
between mergers and RLAGNs is robustly established at 
all redshifts. While both recent simulations and observa¬ 
tions show t hat the galaxies merger rate increases with 
redshift (e . g. iHopkins et ^12006: iGonselice et al.l 120031 : 


iLotz et cd1120111 : iRodriguez-Gomez et al.ll20I5ll . our re¬ 

sults show that there is no evidence that the merger 
fraction for RLAGN is higher at higher z. Of course, this 
should be confirmed through the analysis of larger sam¬ 
ples that may better constrain the merger fractions and 
highlight any possible trends with redshift. But even if 
the uncertainties remain large with the present samples, 
the existence of such a tight relationship between the 
two phenomena means that somehow a RLAGN needs a 
merger (at all redshifts) in order to manifest itself. The 
predominance of major mergers among our radio loud 
samples (at least at z > 1) that is apparent from vi¬ 
sual inspection needs to be confirmed through a more 
quantitative analysis. However, this may imply that one 
of the conditions that must be met in order to trigger 
a RLAGN is that the merger needs to be between two 
galaxies (and thus between two BHs) of similar mass. 
It is well known that radio galaxies are ubiquitously as¬ 
sociated with very massive hosts (> 10^°“^^ Mq, e.g. 
Best et ^1200511 and high-mys SMBH (> 10® M^ ILaor 


20001: Chiaberge fc Marconil 120111: iDuniop et al.i 12005 
Best fc Heckmanll2012f) . Therefore. we expect those ma¬ 


jor mergers to involve high mass objects only. 

6.3. How do mergers affect the central supermassive 
black holes? 


Understanding the details of this issue is central to 
our future work on the subject, but it is beyond the 
scope of this paper. However, we can build on our 
results and speculate on a possible scenario. One of 
the effects of mergers is to lower the specihc angu¬ 
lar momentum of the gas in the gala xy, and thus to 
drive the gas towards the center (e.g. iHernauisH 119891 : 
iHopkins fc QuataertI l201lh . While this may naturally 
happen in gas rich mergers, as those observed in our high- 
z 3CR sample (e.g. iBarthel et al .1120 121: iTadhunter et al.l 
12014 iPodigachoski et al.ll2015ll . the hosts of low redshift 
radio galaxies (in particular those of low radio powers) 
are often relatively gas-poor systems. Therefore, this ef¬ 
fect may play a role at higher redshifts, but it is very 
unlikely to be the ultimate cause of radio loud activity 
in general. Furthermore, tidal effects happen in all merg¬ 
ers, thus merger events should affect both RQ and RL 
AGNs at the same level, contrary to our results. 

Another effect of mergers is to alter the spin and the 
mass of the central black hole. That can be achieved in 
two ways, i.e. either via ac cretion, or via BH-BH merger 
(e.g. iVolonteri et al.l 120131) . The former implies that a 
large amount of gas is driven toward the central region 
of the galaxy, for a significant amount of time. If the 
accretion of matter is coherent, i.e. if the flow of matter 
occurs at a fixed angular momentum axis, that ultimate 
leads to spinning-up the black hole. On the other hand, if 
accretion results from multiple merger events that drive 
the matter towards the black hole from different direc¬ 
tions, the BH is spinned-down. In any case, as already 
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pointed out above, the expected amount of gas in major 
mergers (at least at low redshifts) is probably too small 
to alter the BH spin significantly in these objects. 

In the case of BH-BH merger, the two BH coalesce 
and the resulting object has increased mass and, in most 
cases, a lower spin per unit mass. But there are scenar¬ 
ios in which those events lead to the opposite result. For 
example, it has been shown that a single major BH-BH 
merger, where the ratio between the masses of the two 
involved BHs approaches unity, may generate a spinning 
BH, even if the two merging black h o le are initially not 
spinning i (iHughes fc Blandfordl 120031: [B aker et al.l |2006|: 


iLi et alji201Cll : lGiacomazzo et al.ll201^ SchnittmanI 2013 . 

for a recent review). However, current simulations are 
unable to reproduce BH with dimensionless spin param¬ 
eter gj reater than ^ 0.94 as a result of BH-BH mergers 
alone (jHemberger et al.ll2013fl . 


6.4. The role of rapidly spinning blaek holes 

According; to the so-called spin par adigm 
(jWilson fc ColbertI 119951: iBlandford et al.l 119901 1 , ra¬ 
dio loud AGNs are associated with rapidly spinning 
black holes, while BHs in radio quiet AGN are expected 
to spin less rapidly. Those major mergers that result 
in rapidly spinning black holes may provide the link 
between our observations and the physics behind the 
RLAGN phenomenon as a whole. Clearly, the objects 
that we are seeing in a pre-merger phase, or the very few 
ones that are associated with a minor merger, do not 
fit the above scheme. In those cases, we must assume 
that another major merger event happened in the recent 
past. This does not seem unreasonable, since all of these 
objects lie in over-dense environments, but it should 
be proven by hnding the signatures of that previous 
merger. 

In a framework in which BHs are spun up by ma¬ 
jor BH-BH mergers, we expect a range of result¬ 
ing spin values. Therefore, it is possible that dif¬ 
ferent radio morphologies (and radio powers) are as¬ 
sociated with different BH spin levels. For example, 
only the BHs that spin more rapidly might be able 
to produce the most relativistic jets, in a framework 
in which the jet is powered by energiv extracted from 
the rotating; black hol e fe.g. IBland ford fc ZnaiekI Il977t 
iMcKinnev et al.l l2012t iGhisellini et al.l 120141 1. Such a 
scenario has been explored for X-ray binaries jets, so 
far leading to contrasting results (e . g. [Gardner fc Doni 
2014 iNaravan fc McClintocld 12014 iFender et al.l l20Tu 


Russel^,et,alJ_ 201311. It is interesting to note that 
Fanidakis et al.l ()2nil[ l were able to reproduce the RL 


and RQ AGN populations in the context of galaxy evo¬ 
lution. Their model is based on a scheme that includes a 
bimodal BH spin distribution mainly caused by the the 
combination of accretion of matter onto the black hole 
and by the different type of mergers that galaxies of dif¬ 
ferent stellar mass undergo. 

As pointed out in the previous section, black holes 
could also be spun-up by merger-triggered accretion. 
Such a mechanism could in principle account for our ob¬ 
servations of mergers at z > 1, but only if those are gas- 
rich mergers. However, it would not completely explain 
the origin of radio-loud activity in low redshift, gas-poor 
galaxies (see also the discussion in the next Section). 


6.5. Low redshift radio loud AGNs 

A scenario in which BH-BH mergers are directly 
implied in triggering radio loud AGN activity seems 
to be supported by th e observed properties of low - 
redshift RLAGN host s. iGaoetti fc Balinaverdel (|2006ll : 
Ide Ruiter et al.l (1200511 showed that all radio loud AGNs 
in their samples are hosted by core-galaxies, i.e. galaxies 
that show a flat radial brightness profile in the central 
regions. While this analysis is clearly limited to very 
low redshifts {z < 0.1), where the core radius can be 
resolved in HST images, and to objects that are not 
affected by significant amounts of dust in the central 
kiloparsecs, it very clearly shows that there is a strong 
connection between the presence of a core profile and 
RL activity. Interestingly, one of the most likely expla¬ 
nations for the presence of core profiles is related to a 
major BH-BH mergers, in which the binary BH formed 
during the merger ejects stars fr om the centra l regions 
befor e the two BHs coalesce (e.g. lGrahamll2004 iMerritd 
l2006fl . While the direct physical connection between such 
a phenomenon and radio loud activity is stil l a matter of 
debate (see e.g. iChiaberge fc Marconill2nilh . it provides 
one more piece of evidence that major mergers and radio 
loud AGNs are strictly connected. 

7. CONCLUSIONS 

We derived the merger fraction in samples of z > 1 
Type 2 RLAGN, Type 2 RQAGN, and non- active galax¬ 
ies. We establish that the RLAGN are unambiguously 
associated with mergers (92%('(®J^), while only 38%ti5 
of the RQAGN show evidence for a merger. Non¬ 
active galaxies are statistically indistinguishable from the 
RQAGNs. The comparison with lower redshift samples 
shows that there is no evidence for the fraction of mergers 
in RLAGNs to be dependent on either redshift or AGN 
luminosity. 

Mergers are directly involved in triggering radio-loud 
activity at all redshifts. We speculate that the galaxy 
mergers we observe at z > 1 are responsible for spinning 
up central black holes possibly through mergers of high- 
mass BHs. 

It will be extremely important to study the BH mass 
distribution in samples of merging galaxies of different 
type and to determine whether the BH mass and the 
merger types are related to different type of activity. 

Not all galaxy mergers in our RL samples appear the 
same to visual inspection. It would be important to 
firmly assess the ratio between major and minor mergers 
at different redshifts, and radio power. The images in 
this work show that a larger fraction of objects are in 
a phase of ongoing- or post-major merger in higher red¬ 
shift samples than in t he lower redshift counterparts (see 
also lFlovd et al.l[2nn8[ l. The different HST cameras used 
not only sample different rest-frame wavelengths (near- 
IR at low-z and optical at z > 0.4), but they also have 
significantly different sensitivities. It is extremely hard 
to detect post-merger signatures such as e.g. faint tidal 
tails in the intermediate redshift samples imaged with 
WFPC2. Those features may lie outside the field of view 
in the NICMOS images at low redshift. The only redshift 
range that is covered with homogeneous observations is 
between z=l and 2.5. Clearly needed is a homogeneous 
data set observed with HST-WFC3 or ACS at all red- 
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shifts. 

Detailed studies of the mergers in our radio-loud sam¬ 
ples using integral field spectroscopy, combined with deep 
high resolution imaging in the rest frame IR , should be 
used to study the kinematics of the mergers and mea¬ 
sure a range of parameters such as dynamical masses 
and angular momenta. The imaging part of the project 
is feasible with WFC3, although at high redshifts the 
rest frame wavelengths sampled by such a camera are 
still within the optical bands. Dust obscuration may re¬ 
duce the accuracy of the stellar mass measurements. 

ALMA observations can trace the molecular gas in¬ 
volved in the mergers. It is important to study a large 
sample of these sources at different redshifts, to test 
whether different amounts of gas are driven towards the 
central supermassive black holes and whether that may 
affect the type of radio loud nuclear activity A compari¬ 
son with the mergers observed in radio-quiet AGNs and 
in ULIRGs will elucidate the black hole feeding mecha¬ 
nism in all of these sources, and relationships with AGN 
activity. 

This work is limited to Type 2 AGNs because in Type 1 
AGNs the nuclear light from the central QSO hampers 
morphological studies of the host galaxy. With a better 
knowledge of the WFC3 PSF, such a study is possible 
starting with the sample of SCR QSOs at z > 1. Better 
re sults can be achieved using coronographic observations 
bv IMartel et al.l (1200311 . While the coronographic mode is 
now available on HST only in STIS, JWST with NIRISS 
will allow us to image QSOs spanning a broader range of 
redshift. We hope to be able to study in detail both the 
str ucture and fueling mechanisms in Type 1 objects (see 
e.g. iFord et ^I20l4l . 

In conclusion, our results clearly establish that 
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